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One-particle interchain hopping in a system of coupled Luttinger liquids is investigated by use of 
exact diagonalizations techniques. Firstly, the two chains problem of spinless fermions is studied in 
order to see the behaviour of the band splitting as a function of the exponent a which characterizes 
the ID Luttinger liquid. Moderate intra-chain interactions can lead to a strong reduction of this 
splitting. The on-set of the confinement within the individual chains (defined by a vanishing split- 
ting) seems to be governed by a. We give numerical evidence that inter-chain coherent hopping 
can be totally suppressed for a ~ 0.4 or even smaller a values. The transverse conductivty is 
shown to exhibit a strong incoherent part. Even when coherent inter-chain hopping is believed to 
occur (at small a values), it is shown that the coherent Drude weight is always significantly smaller 
than the incoherent weight. Implications for the optical experiments in quasi-lD organic or high-T c 
superconductors is outlined. 

PACS numbers: 74.72.-h, 71.27.+a 



Recently, the study of strongly correlated fermions con- 
fined to coupled chains has received a great deal of in- 
terest in particular as a way of studying the dimensional 
cross-over from ID Luttinger-like behaviour to 2D. 

Some time ago, Anderson emphasized the crucial dif- 
ference between in-plane and inter-plane (c-axis) trans- 
port observed in copper oxide superconductors jjj. In- 
deed, experimentally the c-axis transport has an anoma- 
lous behaviour ||] in the sense that the transverse conduc- 
tivity has a completely incoherent frequency dependence, 
there seems to be no sizeable Drude-like term (except in 
the optimally doped systems) and <r(u>) is a very slowly 
increasing function of the frequency. This phenomenon 
has been interpreted as an incoherent hopping or as the 
"confinement" of the electrons inside the weakly coupled 
planes. However, for coupled Fermi liquids (FL), Lan- 
dau theory predicts coherent transverse hopping and no 
anomalous transport. Therefore, it has been suggested 
that anomalous transport in the direction of low con- 
ductivity is precisely the signature that the ground state 
(GS) of the two-dimensional (2D) plane itself is not of 
the usual FL type. Various candidates for this state have 
been proposed, such as the marginal Fermi liquid Q or 
the Luttinger liquid (LL) which is generic in one dimen- 
sion (ID) Q. 

Although non- Fermi liquid (NFL) behaviour is thought 
to occur in the planes of high temperature superconduc- 
tors (HTSC), it has been impossible yet to prove the NFL 
nature in 2D except for unrealistic models. However, as 
stated above, it is well known that the generic features of 
correlated ID electrons are not FL-like but rather those 
of a LL and the precise nature (asymptotic behaviour, 
exponents, etc..) of the system can be easily controlled. 
In addition, quasi-lD systems are realized in nature, and 



the problem of coupled chains is of direct relevance there. 
For instance, in the case of the organic superconductors of 
the (TMTSF) 2 X family fj, also known as the Bechgaard 
salts, the high temperature properties are believed to be 
essentially one-dimensional, while the low-temperature 
behaviour is three-dimensional. This cross-over is pre- 
sumably responsible for the anomalies observed in the 
temperature dependence of several quantities (such as the 
2kp contribution to the relaxation rate 0, the ratio of 
the perpendicular conductivity to the parallel one , the 
plasma edge when the electric field is polarized perpen- 
dicular to the chains [|| , . . . ) as well as for the insulating 
behaviour reported for (TMTSF^CICU in the presence 
of a strong enough magnetic field 0. Although a lot 
of work has already been devoted to that problem [MB , 
several aspects of this cross-over have to be better under- 
stood, in particular those concerned with the transport 
properties perpendicular to the chains. 

Hence, from now on and for sake of simplicity, we shall 
only deal with weakly coupled chains. 

The effect of single-particle transverse hopping has pre- 
viously been studied from a renormalization group point 
of view ||]. The notions of interchain coherence or inco- 
herence have emerged from such an analysis. Let us recall 
here that a LL has a different structure from FL: there 
are no quasi-particle like excitations but instead collec- 
tive modes (charge and spin) with different velocities 
which lead to the so-called spin-charge separation; more- 
over, the density of states n(k) has no step-like structure 
at the Fermi level but instead a power law singularity 
n{k) — n{kp) ~ \k — kp\ a defining the parameter a which 
depends on the intra-chain interaction. It turns out that 
the hopping t± is a relevant perturbation when a < 1 [0 . 

However, it has been argued that relevance in that 
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sense was not necessarily a sufficient condition to cause 
coherent motion between chains. This, e.g., can be seen 
from the following model ||; let a system of two sep- 
arated chains be prepared at time t = with a differ- 
ence of AN particles between the two chains. Then, the 
interchain hopping is turned on and one considers the 
probability of the system of remaining in its initial state, 
P(t). Coherence or incoherence can then be defined as 
the presence or absence of oscillations in Pit). 

In Ref. the authors restricted themselves to the 
short time behaviour of P(t) which, they argued, can 
provide valuable informations on interchain coherence or 
incoherence. They found two regimes which depend en- 
tirely on the value of a: the case a < 1/2 exhibited coher- 
ent motion while a > 1/2 showed no signal of coherence. 
However, as has been shown numerically and argued in 
Ref. [fl~Tf] , the small time behaviour can not always dis- 
tinguish between coherent or incoherent regimes. In ad- 
dition, if the initial state is far from equilibrium, i.e. if 
AN/N is not infinitesimal, then mathematical properties 
of the many-body spectrum related to the integrability 
or non integrability of the model will play a central role. 
Indeed, in this case, coherent behaviour is found jflj to 
be a generic property of integrable models (irrespective 
of the value of a). However, it is not clear whether this 
conclusion remains for low energy excitations. 

Another possible issue which has been raised is 
whether the spin degrees of freedom are essential in An- 
derson's confinement mechanism or not JlOj ] . The two 
chain problem with no anomalous exponent (a = 0) but 
with spin-charge separation features was solved exactly 
in |l2| and showed no confinement within each chain i.e. 
an energy separation was found between the bonding and 
anti-bonding states. This result a priori does not con- 
tradict Anderson's conjecture which claims that the be- 
haviour only depends on a. However, it is still not clear 
whether confinement can be found for spinless models 
with large a's and whether the fact that the low energy 
excitations of the spinless fermion chain are collective 
modes is sufficient to, alone, produce confinement in the 
chains. 

This paper is devoted to the study of various aspects 
of interchain coherence in systems of strongly correlated 
spinless fermions. We shall derive several quantities sen- 
sitive to the coherent/incoherent nature of the hopping 
transverse to the chains from exact diagonalizations of 
small systems by the Lanczos algorithm [juj). After a 
short presentation of the models in Sec. I, a precise char- 
acterization of the single chain system in terms of its LL 
parameters is given in Sec. II for models with screened 
interactions extending in space up to third neighbours. 
In Sec. Ill, we study the simplest case of two coupled 
chains. The influence of the interchain hopping on the 
splitting in energy of the singularities appearing in the 
spectral function for kj_ ~ and n is studied, in a sim- 
ilar way as has been done for particles with spin [p"4[ . 
Using different models with different extensions in space 
of the interaction, we investigate whether the physics de- 



pends only on the LL parameter a or not. In order to 
establish possible connections between the splitting and 
transport properties, we consider in Sec. IV a system of 
three coupled chains. Periodic boundary conditions in 
the transverse direction are used in order to derive the 
optical conductivity in this direction, hence, providing a 
direct (and experimentally observable) test of the con- 
finement of the fermions within the chains. 



I. THE MODEL 

We consider here a model of spinless fermions on a 
lattice formed by m chains of length L with a weak inter- 
chain hopping: 
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where (3 labels the chain (j3 = 1, . . . , m), j is a rung index 
(j = 1, . . . , L), Cj_p is the fcrmionic operator which de- 
stroys one fermion at site j on the chain [3, and V{5) is a 
repulsive interaction between two fermions at a distance 
8 (the lattice spacing has been set to one). For con- 
venience, we choose a repulsive interaction of the form 
V(i) = 2V/(i + 1) for i < io, with, more specifically, 
«o = 1,2,3 which corresponds to an interaction extend- 
ing up to first, second and third nearest neighbours (NN) 
respectively. 

In the x and y-directions, we shall use arbitrary bound- 
ary conditions (BC) by threading the system with a mag- 
netic flux <& x and $ y respectively (except for m — 1 or 2 
chains where open BC are used in the y-direction). This 
is realized by a Peierls substitution that modifies the ki- 
netic term, e.g. a twist along the transverse direction is 
realized by 



(2) 



where $j, is the flux measured in unit of the flux quantum 
$0 = hc/e. Similarly, the hopping term in the x-direction 
contains a phase ^j;^x- 

The motivation to introduce flux is two-fold. Firstly, 
as proposed by Kohn |l5| transport properties of a corre- 
lated system can be directly measured from the response 
of the system to a twist in the boundary condition. Sec- 
ondly, our ultimate goal is to extract quantities in the 
thermodynamic limit from finite size scaling analysis. It 
turns out that a simple way to improve the accuracy for 
a fixed system size is to average over the boundary con- 
ditions e.g. over $ y and/or $ x p^Jl7[ |. 
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II. LL PARAMETERS FOR THE SINGLE CHAIN 



In order to characterize the behaviour of coupled 
chains, it is required to, first, compute the parameters 
of a single chain for the same model. Indeed, one impor- 
tant issue is to study whether interchain transport is a 
universal function of the LL parameters only or whether 
it depends on the details of the model. In the case of NN 
interactions, the hamiltonian (known as the t-V model) 
can be mapped onto a spin chain problem by a Jordan- 
Wigner transformation and this is exactly soluble by the 
Bethe ansatz; thus, a is known for each filling Qj. How- 
ever, for extended interactions in space, a numerical in- 
vestigation is necessary with the help of conformal invari- 
ance identities. It turns out that the exponent a can be 
related to simple physical quantities which can be eas- 
ily extracted from standard exact diagonalization results 
using the Lanczos algorithm. For instance the Drudc 
weight 2ttD (D is the charge stiffness) and the charge 
velocity u p are related by ||18|1 



2irD = 2u p K p , 



(3) 



where K p is a universal parameter which determines the 
long distance behaviour of the correlation functions. The 
quantities u p and D can be directly obtained on finite 
systems. The Drude weight for a single chain is given by 
the Kohn formula: 



2ttD = 



1 d 2 (E /L) 
4tt d-K 



(4) 



The charge velocity can be extracted from the difference 
in energy of the two singlet symmetry sectors k = and 
k = 2n / L (for an even number of fermions). K p is then 
obtained from Eq. (||) and the density of state exponent 
(for spinless fermions) can be calculated as Eq] 
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(5) 



Finite size scaling analysis reveals that the 1 /L 2 law for 
the finite size corrections is very well satisfied for K p and 
a so that an accurate determination of the extrapolated 
values can be calculated. Results for a are plotted in 
Fig. | 
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FIG. 1. Exponent a vs V for n = 1/4 and for NN interac- 
tion (♦) and longer range io — 2 (■), io = 3 (•). 

As is expected, the value of a increases with the 
strength of the interaction and with its range in space. 
However, for a given filling, there is a maximum allowed 
value a max for a ^0|. Above this value, Umklapp pro- 
cesses become relevant, a gap opens up in the single parti- 
cle spectrum and the system undergoes a metal-insulator 
transition. The insulating phase is signalled by both 
an opening of the single particle gap and by a vanish- 
ing Drude weight |l5j| . This region is of no interest to 
us in the present study since no transverse hopping oc- 
curs when t± is smaller than the single particle gap. For 
a given density, a acquires its maximum value on the 
metal-insulator transition line. It is important to notice 
that this maximum value depends only on the density 
and, hence, the commensur ability: if n = p/q 7 it is larger 
for larger q. As an example, for n — 1/2 a max = 1/4 
while for n = 1/4 a value as large as a max = 49/16 ~ 3.06 
can be obtained. This fact has motivated our choice of 
n = 1/4 for the following calculations since it gives a 
large range of a values. Typically, for n — 1/4, a max 
is realized by io = 3 and V ~ 7.5. For a shorter range 
interaction, significantly larger values of V are necessary. 



III. INTERCHAIN COHERENCE 

The simplest approach to investigate interchain coher- 
ence is to consider 2 coupled chains, i.e. a 2 x L ladder. 
We proceed along the lines of Ref. 14 1. In the absence of 
interaction, t± leads to bonding and anti-bonding dis- 
persion bands corresponding to transverse momentum 
k± = or 7r respectively, as seen in figure 2. The splitting 
2t± between these bands can be viewed as the signature 
of a coherent transverse hopping. These bands corre- 
spond to a (^-function singularity in the single particle 
hole (electron) spectral function for k < ki<< (k > kp). 
The hole spectral function is defined by, 
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A h (k,oj) 



Im < {<f>a\c 



■is- H + E a 



(6) 



where k = (k, fex) and a similar definition holds for A e 
by exchanging c£. and c k . 

In the case of interacting particles, this <5-function sin- 
gularity is replaced by a power law singularity and the 
elementary excitations are collective modes. Here, we ad- 
dress the issue of the influence of the hopping t± on this 
singularity, in particular we investigate whether a split- 
ting occurs. The choice of the boundary conditions in 
the x-direction is expected to be important in the scaling 
behaviour of various quantities. Periodic or anti-periodic 
BC lead to closed or open shells as can be seen in figure ||. 
We shall consider these two cases separately. 




FIG. 2. Dispersion relation along the chain direction in the 
closed (left) and open (right) shell configurations. Full (open) 
symbols correspond to occupied (empty) states and fj, is the 
chemical potential. 



A. Splitting: closed shell configuration 

In the closed shell configuration of Fig. (||), we can 
add (respectively remove) one fermion on either of the 
two branches (characterized by the transverse momen- 
tum k±) just above (resp. below) the FL and then cal- 
culate the GS energy of this new system. The difference 
between the two values of the energy gives the splitting 
AE between the two bands. 

For interacting particles, AE can be alternatively (and 
more precisely) defined as the energy separation be- 
tween the two low energy peaks in the spectral function 
A/j(k, a;) and A e (k, ui) for k± = and k± = n. For 
a longitudinal momentum k chosen above or below the 
(non-interacting) Fermi wavevector fop = nn the electron 
(hole) spectral function has to be considered. We have 
performed exact diagonalizations on clusters 2 x 4p with 
p=l,...,5atn=l/4to obtain the splitting for differ- 
ent values of the parameters. Results for a 2 x 16 ladder 
at n=l/4 are shown in Fig. (||) (a) and (b) for two mo- 
menta k < kp and k > kp. The results look similar to the 
non-interacting case although the splitting between the 
k±_ = and k± = it structures is reduced. Let us remark 
that a new structure appears for uj < [i and k > kp. 
It is completely absent in the non-interacting case but 
was predicted by Voit (T^j for large a's by using a low 
energy bosonization scheme. Note, however, that such 



calculations can not resolve the fine structure of the LL 
singularity. 
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FIG. 3. Spectral function vs ui for the case V = 5, io = 2 
with t±_ = 0.1. The regions ui < n and u> > /i correspond 
to Ah(k, ui) and A e (k,ui) respectively. The full (dotted) lines 
correspond to k± = n (k± — 0). (a) (resp. (b)) corresponds 
to a momentum below (resp. above) kp. An artificial width 
e has been added to the frequency, u) — > u> + is. 



A careful study of the exact diagonalizations data re- 
veals that AE, for small t±, behaves like at± + bt\ and 
the coefficient a can be calculated accurately. Finite size 
corrections of the order of 1 /L have been found for a so 
that \im t± -,o AE/2t± can be estimated accurately in the 
thermodynamic limit. In Fig. ^, we plot, as a function 
of a, the average between the splittings above and below 
the FL as well as the extrapolated value. The (normal- 
ized) splitting is greatly reduced by the electronic inter- 



4 



actions (measured by a) but it does not seem to vanish 
for a = 0.5. On the contrary, it is compatible with the 
RG approach predicting that t± becomes irrelevant, and 
hence that \im t± ^o AE/2t± vanishes, only for a > 1. 
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the Fermi momentum which is independent of the sys- 
tem size. We thus define the splitting as the difference 
between the electron and hole excitation energies i.e. 
E (N e + l,k e ) - E {N e ) and E {N e ) - E {N e - l,k h ) 
where Eo(N e ) is the reference energy corresponding to 
the absolute GS of the N e = nN electron system. The 
momenta for the electron and hole excitations are fixed, 
k e = (kp,ir) and kfj = (kp,0). Note that, for a fixed 
value of t±, these elementary excitations are not the low- 
est energy excitations in the thermodynamic limit (if t± 
is relevant) since they have a longitudinal momentum hp 
which is defined with respect to the single chain case (i.e. 
for t± = 0). We then naturally define, 

AE = E Q (N e + l,k e )+E (N e - l,k,0 - 2E (N e ). (7) 

For V=0, this expression exactly gives the splitting 2t± 
for any system size. For t± = but finite interaction 
strength, AE is finite. However, Fig. (^) reveals that it 
scales to zero in the thermodynamic limit as expected. 
Fig. Q also shows that an accurate finite size scaling 
analysis can be performed for finite interaction strength 
and finite tx, assuming 1/L finite size corrections. 



FIG. 4. lim ti _,o AE/2t± (=a) vs a for various system sizes 
and extrapolated values. Closed shell configurations have 
been considered. An average between the splitting above and 
below fc_F is performed. L = 8, L — 12 and L — 16 correspond 
to ■, A respectively. 

However, these results must be considered with care 
for two reasons. First, it is clear that AE alone does not 
completely characterize the transverse hopping dynam- 
ics. As seen in Ref. fOfl , the return probability P(t), for 
example, depends on the whole frequency dependence of 
Ah(k, u) (or A e (k, cu)). Since, for large interactions, spec- 
tral weight from the LL singularity is transfered to higher 
excitation energies, then an incoherent behaviour of P(t) 
can be found even though AE remains finite. Secondly, 
it should be stressed that the limit t± — > has been 
taken first before eventually taking the limit L — > oo. In- 
deed, the linear behavior of the splitting with t± is only 
valid in a narrow region which shrinks to zero as 1/L. In 
other words, when t± ~ ttv p /L (v p is some characteristic 
charge velocity) , in the closed shell configuration, the ex- 
citation energy at k± = n for a longitudinal momemtum 
just below kp crosses the chemical potential so that the 
splitting cannot any longer be defined as some energy dif- 
ference between many-body states corresponding to the 
same number of particles. To avoid such complications, 
we now consider the open shell configuration of Fig. ^j. 



B. Splitting: open shell configuration 

As can be seen in Fig. (||), in the open shell con- 
figuration, one can add or remove a particle exactly at 
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FIG. 5. AE vs 1/L for t± = (open symbols) and t± = 0.3 
(solid symbols). Interactions up to next nearest neighbors 
(io = 2) with different strength V have been considered. 

The extraplolated values (L = oo) of the ratio A/2t± 
are plotted as a function of the exponent a of the io = 2 
model (i.e. with interactions extending up to second 
nearest neighbors) in Fig. (|^) for fixed interchain hop- 
ping amplitudes tj_. A strong reduction of this ratio 
for increasing a indicates that intrachain repulsion has a 
drastic influence in prohibiting interchain coherent hop- 
ping. Our data suggest that there is a critical value 
a* (t± ) such that interchain hopping becomes incoherent 
for a > a*(t±). As expected, a*(t±) increases with tx- 
Even for finite (but small) t± , a* (t± ) can be as small as 
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0.4. The limiting value a*(t± = 0) = ao, although not 
accurately given by our method, is probably significantly 
smaller than 0.4. This is indeed clear from Fig. (Q) where 
the same ratio A/2t± is plotted as a function of t± for 
various interactions (iq = 2 model). For large interac- 
tions like V — 3.25 (a ~ 0.38), there is a critical value 
t*i(a) of t± below which incoherent transverse hopping 
takes place. This corresponds to the case a > uq. For 
smaller interaction like V = 1 (a ~ 0.08) our data are 
consistent with A/2t± approaching a finite value when 
t±_ —> 0. Hence, transverse hopping remains coherent in 
this case which corresponds to a < ao . For intermediate 
interactions like V = 2 (a ~ 0.22), our data are not con- 
clusive. However, very small values of ao like 0.2 (or even 
smaller) are not inconsistent with our numerical analysis. 




into account, the spin-charge separation that occurs in 
ID should suppress even further the coherent transverse 
hopping. At a qualitative level, this can be understood 
from the fact that only real electrons can hop from one 
chain to the next and this is believed to become more 
difficult in the presence of spin-charge separation as has 
been suggested by Anderson Q. Although dealing with 
a different filling n = 1/3 and with particles with spin, 
the same qualitative behaviour is found numerically in 
Rcf. |L4| ; however, the decrease of the splitting with a is 
stronger. 
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FIG. 7. Extrapolated (L = oo) values of AE/2t± vs t± in 
the io = 2 spinless fermion model with several values of V 
(indicated on the figure). 



FIG. 6. Extrapolated (L = oo) values of AE/2t± vs a for 
several values of t± (indicated on the figure) and in the io = 2 
spinless fermion model. 

It is important to notice that the previous extrapola- 
tions carried out for open shell configurations and finite 
t± give very different results than those shown in Fig. (|4|). 
In particular, Fig. (||) suggests that the critical value ao 
is significantly smaller than the value of 1 predicted by 
the RG analysis. One possible explaination is that the 
limits tj_ — > and L — > oo do not commute with each 
other. This scenario is supported by the fact that, in the 
closed shell configuration, the linear dependance of the 
splitting AE with t± is limited, for incresing system size, 
to a narrower and narrower range of order X/L. There- 
fore, an extrapolation at finite t± has to be realized. It is 
interesting to note that the RG approach seems, on the 
contrary, to reproduce the lmiL—, 00 {lmit ± ^o(AE/2t±)} 
data. 

We shall finish this section with a brief comment on 
the role of the spin degrees of freedom which have not 
been considered here. It is clear that, when spin is taken 



IV. TRANSPORT PROPERTIES 

The previous study suggests that the interaction tends 
to confine the electrons within the chains, although no 
complete confinement seems to occur at small a values 
where AE ^ 0. A better understanding of this phe- 
nomenon can be achieved by investigating the transport 
properties along the y-axis (inter chain) and more pre- 
cisely the transverse optical conductivity which is the 
linear response of the system to a spatially uniform, time 
dependent electric field in the transverse direction. For 
such a study, a torus geometry is needed (to > 3) so that 
a current can flow around the loop in the y-direction. 
One of the main advantage of the optical conductivity is 
that it can directly be measured experimentally. 

The real part of the optical conductivity can be written 
as a sum of two parts, 

* vy {u) = 2nD vy S{w) (8) 

where D yy is the charge stiffness in the y-direction which 
defines the Drude weight 2-nD vv . Note that an important 
f-sum rule M can be used to check the numerical results, 
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rT v . 



(9) 



where the expression between brackets is the mean value 
of the transverse kinetic energy in the ground state. 

The actual calculation of the frequency dependence of 
(jy y 9 (uj) will be carried out later on and we first concen- 
trate on D yy and on the total sum rule. As originally 
noted by Kohn D yy , which measures a transport 
quantity, can be obtained from the dependence of the 
ground state energy Eq on Q y as 



2irD m ($ v ) 



m 2 d 2 {E /N) 



4tt 



d$2 



(10) 



where N = mL is the number of sites. This is, in all 
points, very similar to the previous derivation of the 
charge stiffness D for the ID rings. 

The previous quantities have been calculated numeri- 
cally on finite 3 x L lattices using the Lanczos algorithm. 
We have chosen a quarter filled band so that the extrap- 
olation of results for L = 4, 8 and 12 is possible. From 
now on, we shall restrict ourselves to PBC or ABC in 
the x-direction. In most cases, an average over <fr y is re- 
alized jl6 17 in order to mimic the case of many parallel 
chains i.e. we calculate, 



,1/2 

{Dyvhy = / d%D yy (%). 

J -1/2 



(11) 



This is done by calculating the ground state energy for a 
(large) discrete set of flux values (still making use of the 
translational invariance). A typical curve for Eo($ y ) is 
shown in Fig. (|J) and reveals that level crossings occur 
as a function of <5> y . As for the ladder case, the scaling 
behaviour depends crucially on the choice of the bound- 
ary conditions along the x-direction. In the following, we 
shall study two different cases separately. 
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FIG. 8. GS energy vs & y for a 3 x 8 systems with V — 2, 
io = 2, £_l = 0.7 and closed shell configuration at quarter 
filling . The thick curve is the absolute GS energy. 



A. Charge stiffness: closed shell configuration 

In the non-interacting case of a 3 x L system with peri- 
odic BC in both directions, t± leads to three branches in 
fc-space separated by an energy of order ~ t ± (see Fig. ||) . 
Note that for special values of $ y (e.g. $ a = 0, |), two 
of them are degenerate. We first consider the closed shell 
configuration where, for a given momentum k along the 
chain direction (and sufficiently small t±), the three pos- 
sible states with different k± momenta (0, or — ^) are 
either fully occupied or completely empty. In this case, 
there is a critical value of t ± for which, for an optimum 
choice of $ y , one can move a fermion from one branch 
to the next with no en ergy cost. A straightforward cal- 
culation gives ij_ = \/2/3 sbx(ir/L). Below this value, 
in the y-direction the bands are filled and therefore no 
transport can occur in the transverse direction and the 
Drude weight is vanishing. Note that this quantity t*j_ is 
directly related to the splitting of the bands. 




FIG. 9. Dispersion relations as a function of k for closed 
and open shells. 

More generally in the interacting case, the same be- 
haviour is observed as can be seen in Fig. [l^ where we 
plot the Drude weight as a function of t± for a given 
size. This is a typical behaviour from which one can 
define some precise cross-over value t*j_(L) where a tran- 
sition occurs (D yy increases suddenly). Of course, this 
behaviour is a finite size effect and we must be careful to 
extract thermodynamic results properly. 
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FIG. 10. Drude weight in arbitrary units as a function of 
t± for the case V — 2, io = 3 on the 3x4 and 3x8 lattices. An 
average over $ H has been performed according to Eq. (pi]). 

Two scenarios can occur. Firstly, t*j_(oo) might remain 
finite. In this case, no coherent transport occurs between 
the chains for sufficiently small t±. This is the case, for 
example, for very strong interactions in the insulating 
phase. This is also expected in the LL phase for a > a^. 
Indeed, we expect that AE = (see sec. Ill) would 
imply D yy = (for L = oo). Although it is difficult 
to prove this behaviour numerically, it is not incompati- 
ble with our numerical results. The second possibility is 
that t x (L) — > when L — > oo. In this case, we expect 
that the non-interacting picture should be approximately 
valid, at least for not too large an interaction. In other 
words, according to the non- interacting picture which we 
discussed above, we expect t x (L) to be directly related 
to the splitting AE calculated in the previous section. 



Lt 



perp 




FIG. 11. Lt* x [L) as a function of 1/L for the 
non-interacting case (•) and for finite interaction io = 2 and 
V = 1,2, 3.25, 5 (♦) 

We proceed as follows; we assume a connection be- 
tween t* ± (L) and the splitting AE and show that this 
hypothesis is, indeed, consistent with the numerical data. 
Let us suppose that t* ± (L) fulfills, 



2tt 

4> H L 



(12) 



where v p is some charge velocity. Thus, t* ± (L) should 
be of order 1/L, more precisely it should behave like 
t* ± (L) = A/L + B/L 2 . Finite-size scaling can be per- 
formed by plotting Lt* ± (L) vs 1/L. As seen in Fig. (ll|), 
this scaling law seems to be well satisfied for small in- 
teraction strength so that, in this case, a finite-size ex- 
trapolation of A is possible. On the other hand, a di- 
verging value of Lt± (or A) would mean that is 
in fact finite and, hence, completely incoherent hopping 
occurs. This is likely to occur for the largest values of 
V we have considered, such as V = 3.25 and V = 5 and 
it can not be completely excluded for smaller values of 
V. According to the qualitative argument of Eq. (p"2]), 
if AE(t x ,^ y ) is linear in t± for small t±, then 1/A is 
expected to be directly proportional to the extrapolated 
value of AE. In order to check this point, it is conve- 
nient to normalize A with respect to the non-interacting 
case, A — ^/2/37r for n = 1/4. The quantity A/Aq is 
plotted as a function of the Luttinger parameter a for dif- 
ferent values of the interaction in Fig. [l2| and compared 
to the extrapolated value of AE obtained for closed shell 
configurations. Firstly, we remark that, although some 
error bars are large, the behaviour of t* ± seems to be 
quite similar for different kinds of interactions. In other 
words, a seems to be the only parameter controlling the 
behaviour of 1/A. Secondly, we observe some discrepan- 
cies between 1/A and AE. They can be attributed (i) 
to the crudeness of the picture we have developed only 
at a qualitative level, (ii) to the strong dependence of 
the charge velocity v p in Eq. ( [l2] ) with the interaction 
and (iii) more importantly, to the fact that the previ- 
ous estimation of the splitting realized in the closed shell 
configuration case seems not to be relevant for finite t±. 
Indeed, the second estimation of Sec. III.B of the ratio 
AE/2t±, realized for open shell configurations and finite 
t_i_, gives, for small a, smaller values in better agreement 
with Ao/A. Moreover, the above extrapolation of t* ± (L) 
assumes t*j_(L = oo) = 0, i.e. coherent interchain trans- 
port. However, as mentioned above, the possibility of a 
small finite value of t x (L = oo) cannot certainly be ex- 
cluded when a > 0.2 which would be in complete agree- 
ment with the results for AE/2t± obtained in the open 
shell case. 
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FIG. 12. Extrapolated values of Ao/A as a function of a 
for different range of interactions io = 1 (♦), io = 2 (■) and 
io = 3 (•). For comparison, the extrapolated values of AE of 
Fig. (Q) are also reproduced here (open circles). 

At this point, we are still left with an incomplete 
picture. For not too large an interaction (a < oto), 
our results suggest that Dyy is finite for small trans- 
verse hopping. Fig. ( |Io| ) suggests that the slope 
d(D vy )/dtx \t±=t* seems to vanish when L — > oo. There- 
fore, a tj_ behavior of the Drude weight with t± is cer- 
tainly possible as in the non-interacting case. To have 
a better understanding of this behaviour we now use a 
different choice of the boundary conditions along x. 

B. Charge stiffness: open shell configuration 

By choosing adequate boundary conditions along x 
(<5 X = or & x = ir depending on the length L), open 
shells can be realized as seen in Fig. (^). In this case, 
even for a finite system, the Drude weight and the to- 
tal kinetic energy remains finite down to vanishing ij_ 
as in the non-interacting case. We first consider a fixed 
value of the flux Q y . It is in fact convenient to choose 
<fr y corresponding to the lowest GS energy, i.e. such that 
= 0, since this value (in fact $ y ~ 0.25) is almost in- 
dependent of the interaction and of the system size. The 
corresponding Drude weight is shown in Fig. [l3]for 3x4, 
3x8 and 3 x 12 systems as a function of t±. Finite size 
effects are found to be already weak for the two largest 
cases L = 8 and L — 12. The Drude weight is clearly 
strongly suppressed compared to the V = case |^2| . 



FIG. 13. irDyy for a fixed value of <& y (see text), V — 2, 
io = 3 and different system sizes L — 4,8, 12 from up to down 
(open symbols) compared to the averaged data of L — 8 (filled 
squares). The non- interacting case (for a 3 x 8 torus) is shown 
as a dotted line. 

We now consider an average over <E> y and we shall de- 
note by I yy = (J °° <jy y (u))duS)<s> y the total sum which is 
obtained by computing the transverse kinetic energy and 
averaging over <E> y . As seen in Fig. [l^, the values obtained 
for the Drude weight by averaging the 3x8 data over Q y 
are very close to the ones obtained on the 3 x 12 cluster 
at constant flux. Therefore, it is technically avantageous, 
as far as CPU time is concerned, to consider smaller sys- 
tems and perform a flux average. Note that, however, if 
the number of coupled chains is kept fixed (here m = 3), 
even in the limit L — > oo, we expect (D yy )q> y ^ D yy (<& y ) 
(since the thermodynamic limit is not taken in the direc- 
tion of the flux) . Qualitatively, the averaging procedure 
mimics many coupled chains. In any case, when con- 
finement within each of the individual chains starts to 
occur we do not expect crutial differences between the 
cases of 3 or of an infinite number of coupled chains (if 
only the t± term couples the chains). In Fig. [l4| we plot 
the Drude weight and the total sum rule as a function 
of ij_ for a moderate electronic interaction. We observe 
that the behavior of these two quantities is not incom- 
patible with the t\ law of the non- interacting case [ p2[ . 
However, the intrachain interaction has drastic effects. 
Firstly, the total sum rule is strongly reduced compared 
to V = (shown as a reference on Fig. HI). Secondly, 
it is found that n{D y y)$ and I yy behave differently. In 
other words, for small t± the main fraction of the weight 
lies in the incoherent part. 
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FIG. 14. it{Dyy)<s> y (■) and I yy (•) as a function of t± 
for a 3 x 8 system and V — 2, io = 3. For comparison, the 
V — case is also displayed as a dotted line. 



FIG. 15. Ratio of the Drude weight to the total conductiv- 
ity as a function of t± for NN and longer range interaction 
with the same NN magnitudes V = 2. 



To be more quantitative, it is convenient to consider 
the ratio 

r ADyyh y mVggA ,,_ T > (13) 

which, as can be seen from the sum rule Eq. (||), cor- 
responds to the relative part of the Drude weight in the 
total optical conductivity. We have plotted the ratio r for 
the same 3x8 torus and for the three interaction ranges, 
?'o = 1, 2 and 3. The common important feature of these 
data is that the ratio r decreases as t± goes to and can 
become rather small, typically smaller than 20%. Unfor- 
tunately, we think our data become somewhat unreliable 
for very small t± (let's say t± < 0.1) so that the behavior 
of the ratio r when t± — > cannot be accurately deter- 
mined. However, our numerical estimates should give the 
correct trend in the range 0.1 < t± < 0.4. 



These results should have very important consequences 
on the experimental side. Indeed, our results predict, 
for small t±, anomalous transport perpendicular to the 
chains even when D yy does not completely vanish (if 
a < cto). Spectral weight is suppressed predominantly 
from the coherent part of the conductivity. In order to 
confirm this behaviour, the frequency dependent optical 
conductivity has been calculated directly and results are 
discussed in the following. 

C. Optical conductivity 

The frequency dependence of the transverse optical 
conductivity can be calculated by use of the Kubo for- 
mula, 

»W = ^E ^f^f 6{u - (E n - Eo)) (14) 

where j y is the transverse current operator and the sum 
runs over all the excited states. So far, we restrict our- 
selves to the open shell configuration used in the preced- 
ing Section. o r yy 9 {u) can be calculated exactly on the 
same finite clusters by a continuous fraction expansion 
generated by use of the Lanczos method. Firstly, for a 
given value of <E> y , one computes the ground state. Here, 
we choose the absolute GS which carries no current in 
the y-direction. Then, by applying the transverse cur- 
rent operator on this state, one generates a new vector 
which serves as the starting point of another Lanczos pro- 
cedure. Eventually, the tridiagonal form of the hamilto- 
nian in this new basis is used to compute the continued 
fraction expansion of the regular part of a yy (uj). 
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FIG. 16. Transverse optical conductivity vs frequency for 
a 3 x 8 system at quarter filling with V — 2, io = 3 and 
i_L = 0.1. The Drude 5-function has been represented with 
the same small imaginary part e = 0.04. 

In the free case, the current operator commutes with 
the Hamiltonian and therefore, the conductivity only 
contains a Drude peak. However, as can be seen on 
Fig. [l6], for finite interaction strengths and in the 3x8 
cluster, a pronounced structure appears at finite fre- 
quency for any small value of t±. In order to determine 
more accurately the position in energy of the weight, we 
have computed the first moment of the distribution, 

/>oo />oo 

( W ) = / ff^Hwdw/ / &™{w)du (15) 
Jo+ Jo+ 

which is expected to behave smoothly with the various 
parameters. We observe that it has a finite limit as t± 
goes to (see Fig. [i"7| ). This means that, once t± is turned 
on between the chains, weight immediately appears pre- 
dominantly at finite frequencies. This typical frequency 
increases with the strength and with the range of the 
interaction. This is clearly a signature of some form of 
incoherent perpendicular transport. 



FIG. 17. First moment of the conductivity as a function of 
t± for V = 2 and different interaction ranges. 

V. CONCLUSIONS 

In this work, different approaches to interchain coher- 
ence have been investigated. As a first step, we have 
focussed on the energy splitting generated by the trans- 
verse hopping in the dispersion relation of the LL collec- 
tive modes. By finite size scaling analysis, we have shown 
that this splitting was monitored by the LL parameter 
a. However, incoherent interchain hopping is found for 
much smaller values of a than those predicted by the RG 
calculations || . In the second part of this work, we have 
attempted to make the connection between the previous 
energy splitting and transverse transport properties. In 
the regime where t± is still relevant (a < ao), the most 
important results are that (i) the Drude weight and the 
total optical sum rule grow less rapidly with t± than in 
the non interacting case and (ii) the Drude weight be- 
comes significantly smaller than the total sum rule when 
the intrachain interaction is turned on. Hence, even when 
the Drude weight remains finite (when t± is relevant), 
transverse transport is predominantly incoherent in the 
small t± regime. How small t± needs to be so that this 
regime is observed depends on the strength of the in- 
teraction. Typically, for a ~ 0.2, strong suppression of 
coherent transport occurs up to t±_ ~ 0.15. This phe- 
nomenon could explain the anomalous transport which 
is observed experimentally. 
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